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In planar nano- magnetic devices magnetization direction is kept close to a given plane by the large 
easy-plane magnetic anisotropy, for example by the shape anisotropy in a thin film. In this case 
magnetization shows effectively in-plane dynamics with only one angle required for its description. 
Moreover, the motion can become overdamped even for small values of Gilbert damping. We 
derive the equations of effective in-plane dynamics in the presence of spin-transfer torques. The 
simplifications achieved in the overdamped regime allow to study systems with several dynamic 
magnetic pieces ("free layers"). A transition from a spin-transfer device with a static polarizer to 
a device with two equivalent magnets is observed. When the size difference between the magnets 
is less than critical, the device does not exhibit switching, but goes directly into the "windmill" 
precession state. 
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I. INTRODUCTION 

The predictior^ii^ and first experimental 
observations^iii^iSiili^ of spin-transfer torques opened a 
new field in magnetism which studies non-equifibrium 
magnetic interactions induced by electric current. Since 
such interactions are relatively significant only in very 
small structures, the topic is a part of nano-magnetism. 
The current-induced switching of magnetic devices 
achieved through spin-transfer torques is a candidate 
for being used as a writing process in magnetic random 
access memory (MRAM) devices. The MRAM memory 
cell is a typical example of a spintronic device in which 
the electron spin is used to achieve useful logic, memory 
or other operations normally performed by electronic 
circuits. 

To produce the spin-transfer torques, electric currents 
have to flow through the spatially non-uniform mag- 
netic configurations in which the variation of magneti- 
zation can be either continuous or abrupt. The first case 
is usually experimentally realized in magnetic domain 
walls4i^i^iii Here we will be focusing on the second case 
realized in the artificially grown nano-structures. Such 
spin-transfer devices contain several magnetic pieces sep- 
arated by non-magnetic metal spacers allowing for arbi- 
trary angles between the magnetic moments of the pieces. 
Magnetization my vary within each piece as well, but that 
variation is usually much smaller and vanishes as the size 
of piece is reduced, or for larger values of spin-stiffness 
of magnetic material. The typical examples of a system 
with discrete variation of magnetization are the "nano- 
pillar" devices^ (Fig.lA). Their behavior can be reason- 
ably well approximated by assuming that magnetic pieces 
are mono-domain, each described by a single magnetiza- 
tion vector M{t) = Msn{t) where ft is the unit vector 
and Ms is the saturation magnetization. The evolution 
of n{t) is governed by the Landau-Lifshitz-Gilbert (LLG) 
equation with spin-transfer termsi^ii^ 

It is often the case that magnetic pieces in a spin- 



transfer device have a strong easy-plane anisotropy. For 
example, in nano-pillars both the polarizer and the free 
magnetic layer are disks with the diameter much larger 
than the thickness. Consequently, the shape anisotropy 
makes the plane of the disk an easy magnetic plane. In 
the planar devices'^^ built from thin film layers (Fig. IB) 
the shape anisotropy produces the same effect. When 
the easy-plane anisotropy energy is much larger then all 
other energies, the deviations of n{t) from the in-plane di- 
rection are very small. An approximation based on such 
smallness is possible and provides an effective description 
of the magnetic dynamics in terms of the direction of the 
projection of n{t) on the easy plane, i.e. in terms of one 
aziniuthal angle. In this paper we derive the equations for 
effective in-plane motion in the presence of spin-transfer 
effect and discuss their use by considering several exam- 
ples. 

In the absence of spin-transfer effects the large easy- 
plane anisotropy creates a regime of overdamped mo- 
tion even for the small values of Gilbert damping con- 
stant a <^ iJ^ In that regime the equations simplify 
further. Here the overdamped regime is discussed in 
the presence of electric current. The reduction of the 
number of equations allows for a simple consideration 
of a spin-transfer device with two dynamic magnetic 
pieces. We show how an asymmetry in the sizes of these 
pieces creates a transition between the polarizer-analyzer 
("fixed layer - free layer") operation regim o^i^i^^i^^ and 
the regime of nearly identical pieces where current leads 
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FIG. 1: Planar spin-transfer devices 
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not to switching, but directly to the Slonczewski "wind- 
mill" dynamic statci^ Finally, we point out the limita- 
tions of the overdamped approximation in the presence 
of the spin-transfer torques. 



II. DYNAMIC EQUATIONS IN THE LIMIT OF 
A LARGE EASY-PLANE ANISOTROPY 



Magnetization dynamics in the presence of electric cur- 
rent is governed by the LLG equation with the spin- 
transfer term.~ii2, For each of the magnets in the device 
shown on Fig. lA 
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where s{t) is the unit vector along the instantaneous 
magnetization of the other magnet and the spin-transfer 
magnitude 
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is proportional to the electric current /. Here e is the 
(negative) electron charge, so u is positive when electrons 
flow into the magnet. Due to the inverse proportionality 
to the volume V, the larger magnets become less sensi- 
tive to the current and can serve as spin-polarizers with 
a fixed magnetization direction. As for the other pa- 
rameters, 7 is the gyromagnetic ratio, g{P, {ft ■ s)) is the 
Slonczewski spin polarization factor— which depends on 
many system parameters, ^^'^^ and a is the Gilbert damp- 
ing which also depends on n and s when spin pumpingi^ 
is taken into account. We will restrict our treatment to 
the constant g and a to focus on the effects specific to 
the strong easy plane anisotropy. 

In terms of the polar angles {6, (f)) the LLG equation 
([1]) has the form 
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where the tangent unit vectors eg and are defined in 
Appendix [X] 

We will consider a model for which the energy of a 
magnet is given by 
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with K± being the easy-plane constant, Er being the 
"residual" in-plane anisotropy energy and z-axis directed 
perpendicular to the easy plane. The limit of a strong 
easy-plane anisotropy is achieved when the maximal vari- 
ation of the residual energy is small compared to the 
easy-plane energy, Ai?,. -C K±. In this case = tt/2 + 69 
with (56* < 1. 



To estimate SO, consider the motion of magnetization 
initially lying in-plane off the minimum of Er and neglect 
for the moment the spin-transfer terms in Eq.Q. Mag- 
netization starts moving and a certain deviation from the 
easy plane is developed. For the estimate, assume that 
the energy is conserved during this motion (the presence 
of damping will only decrease SO). Then 
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We can now linearize the right hand sides of equations Q 
in small 56. On top of that, some terms on the left hand 
sides of @ turn out to be small and can be discarded. 
Indeed, taking into account the smallness of a one gets 
the estimates 



7 dEr 
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Consequently ^ (py^ AEr/K± ^ (j) and (j> ^ aO, there- 
fore the second term on the left hand side of the second 
equation of the system ([3]) can be discarded. No simpli- 
fication happens on the left hand side of the first equa- 
tion, where 9 and acj) can be of the same order when 



< 



^AEr/Kj 



Putting the spin-transfer terms back we get the form 
of equations in the limit of large easy-plane anisotropy: 
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Expressions for the scalar products in ([6]) in terms of 
polar angles are given in Appendix [X] 

The second equation shows that 59 can be expressed 
through ((/),</)). Small out-of-plane deviation becomes a 
"slave" of the in-plane motion;^ We get 
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The term with the second time derivative </) decreases 
with increasing Ki_ . As pointed out in Ref. [3, in the 
absence of spin-transfer this term can be neglected when 
K± > AEr/a^. Mathematically this corresponds to a 
transition from an underdamped to an overdamped be- 
havior of an oscillator as the oscillator mass decreases. 

With spin-transfer terms the overdamped approxima- 
tion gives an equation 
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where ^ = uMg/ {-yK^}. The range of this equation's 
validity will be discussed in Sec. II VI The scalar products 
in Eq. ([5]) have to be expressed through the polar angles 
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{ds{t), (t>s{t)) of vector s, and linearized with respect to 69 
(see Appendix, Eq. IA4p . which is then substituted from 
Eq. ([6]). Finally, the equation is linearized with respect 
to small spin-transfer magnitude u. We get: 
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describing the in-plane overdamped motion of an ana- 
lyzer with a polarizer pointed in the arbitrary direction. 
Next, we show how some known results on spin-transfer 
systems are recovered in the approximation 

Consider the device shown on Fig. lA and assume that 
the first magnet is very large. As explained above, this 
magnet is not affected by the current and serves as a fixed 
source of spin-polarized electrons for the second magnet 
called the analyzer, or the "free" layer. The magneti- 
zation dynamics of the analyzer is described by Eq. ([3]). 
The case of static polarizer is extensively studied in the 
literature. 

First, consider the case of collinear switching, exper- 
imentally realized in a nano-pillar device with the ana- 
lyzer's and polarizer's easy axes along the x direction: 
Er = {l/2)K\ism'^(j), s = (1,0,0).^ Using Eq. ® with 
9s — n/2, (j)s = we get 
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Without the current, there are four possible equilibria 
of the analyzer. Two stable equilibria are the parallel 
((/) = 0) and anti-parallel (0 — tt) states. Two perpen- 
dicular equilibria (0 — ±7r/2) are unstable. Lineariz- 
ing Eq. (jlOp near equilibria one finds solutions the form 
5(j){t) exp^Lut) with eigcnfrcquencies 
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The equilibria are stable for a; < and unstable other- 
wise. Thus the parallel state is stable for ^ > —a/2, the 
antiparallel state is stable for ^ < a/2, and the perpen- 
dicular states cannot be stabilized by the current. These 
conclusions agree with the results of Refs. iQIll The 
stability regions are shown in Fig. 

Note how Eq. (fTO|) emphasizes the fact that spin- 
transfer torque destabilizes the equilibria by making the 
effective damping constant ae// = a + 2£^ coscf) negative, 
while the equilibrium points remain a minimum of the en- 
ergy Er ■ Any appreciable influence of the current on the 
position and nature (minimum or maximum) of the equi- 
librium can only be observed at the current magnitudes 
1/a times larger than the actual switching currentJ^ 
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FIG. 2: Stability regions for systems with static (A) and 
dynamic (B) polarizers as a function of applied current, 
C = g{P){h/2VK±)I/e oc /. 



Second, consider the case of magnetic fan.^ Here the 
easy axis of the polarizer is again directed along x, but 
the polarizer is perpendicular to the easy plane: s — 
(0,0,1), 9s — 0. This arrangement is known to produce 
a constant precession of vector n. Eq. ^ gets a form: 
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for |7i| < jK\\/ {2Ms) the current deflects the analyzer 
direction from the easy axis direction. For larger values 
of u there is no time-independent solution. The angles 
(j) grows with time which corresponds to ft making full 
rotations. At |u| ^ 'yK^/{2Ms) the rotation frequency 
of the magnetic fan is given by w ~ u/a. 



III. DEVICE WITH TWO DYNAMIC 
MAGNETS (TWO "FREE LAYERS") 

No let us assume that both magnets in Fig. lA have 
finite size. Each magnet serves as a polarizer for the other 
one. Without approximations, the evolution of two sets 
of polar angles {9i, i = 1,2 is described by two LLG 
systems of equations 
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where j means the index not equal to i and no summation 
is implied. 

We now apply the overdamped, large easy-plane 
anisotropy approximation to both magnets. Equation ^ 
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for each magnet is further simplified since for the magnet 



i the angle 9^ — 6j — 7r/2 



<^ 1. Expanding 



in small 69 j and using the slave condition ([6]) for 69 j 
with (s • e<j4) = (n'--'^ • el*'') expanded in both small angles 



(see Eq. ()A5P ) we get the system: 
(ctj + 2^ji cos((/)j - 
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with = UjiMsi/ {'~fK±). It was assumed that K± is 
the same for both magnets. 

The spin-transfer torque parameters U21 and U12 have 
opposite signs and their absolute values are different due 
to different volumes of the magnets, according to Eq. ([2]). 
We assume Vi > V2 and denote U12 = u, U21 = —eu. 
The larger magnet experiences a relatively smaller spin 
transfer effect, and the asymmetry parameter satisfies 
< e < 1. In general, material parameters ai^2, Afsi,2 
and magnetic anisotropy energies E^^''^"> of the two mag- 
nets are also different, but here we focus solely on the 
asymmetry in spin-transfer parameters. Both E^^"> and 
are assumed to be given by formula (|4|) with the 
same direction of in-plane easy axis. The situation can 
be viewed as a collinear switching setup with dynamic 
polarizer. Equations (|13p specialize to 
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Next, we study the stability of all equilibrium configu- 
rations ((/>i,02) of two magnets. There are four equilib- 
rium states that are stable without the current: two par- 
allel states along the easy axis (0,0) and (tt, tt), two an- 
tiparallel states along the easy axis (0, tt) and (tt, 0). Four 
more equilibrium states have magnetization perpendicu- 
lar to the easy axis and are unstable without the current: 
(±7r/2, ±7r/2). Once again, since spin-transfer does not 
depend on the relative direction of current and magneti- 
zation, the configurations which can be transformed into 
each other by a rotation of the magnetic space as a whole 
behave identically. Thus it is enough to consider four con- 
figurations: (0,0), (0,7r), (7r/2,7r/2), and (7r/2, -7r/2). 
We linearize equations (|14p near each equilibrium and 
search for the solution in the form 6(f)i ^ expicut). The 
eigenfrequencies are found to be: 
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The state is stable when both eigenfrequencies are 
negative. We conclude that initially unstable states 
(7r/2,±7r/2) are never stabilized by the current, while 
the (0, 0) and (0, tt) state remain stable for 

These regions of stability are shown schematically in 
Fig. [2J3 in comparison with the case of static magnetic 
polarizer (Fig. Wh) which is recovered at e — > 0. 

As the size of the polarizer is reduced, the asymme- 
try parameter e grows. The stability region of the an- 
tiparallel state acquires a lower boundary ^ = ^a/{2e). 
Up to e = 1/2, this boundary is still below the lower 
boundary of the parallel configuration stability region. 
Consequently, the parallel configuration is switched to 
the antiparallel at a negative current ^ = — q;/(2(1 — e)). 
The system then remains in the antiparallel state down 
to ^ = —a/{2e). Below that threshold no stable configu- 
rations exist, and the system goes into some type of pre- 
cession state. This dynamic state is related to the "wind- 
mill" state predicted in Ref. ^ for two identical magnets 
in the absence of anisotropics. Obviously, here it is mod- 
ified by the strong easy-plane anisotropy. 

The e = 1/2 value represents a transition point in the 
behavior of the system. For 1/2 < e < 1, the stability 
region of the parallel configuration completely covers the 
one of the antiparallel state. A transition without hys- 
teresis now happens at ^ = — a/(2(l — e)) between the 
parallel state and the precession state. If the system is 
initially in the antiparallel state, it switches to the par- 
allel state either at a negative current ^ — — q;/(2(1 — e)) 
or at a positive current ^ = a/ 2, and never returns to 
the antiparallel state after that. 



IV. CONCLUDING REMARKS 

We studied the behavior of planar spin-transfer devices 
with magnetic energy dominated by the large easy-plane 
anisotropy. The overdamped approximation in the pres- 
ence of current-induced torque was derived and checked 
against the cases already discussed in the literature. In 
the new "dynamic polarizer" case, we found a transition 
between two regimes with different switching sequences. 
The large asymmetry regime is similar to the case of 
static polarizer and shows hysteretic switching between 
the parallel and antiparallel configurations, while in the 
small asymmetry regime the magnets do not switch, but 
go directly into the "windmill" precession state. 

We saw that the current-induced switching occurs 
when the effective damping constant vanishes near a par- 
ticular equilibrium. This makes the overdamped approxi- 
mation inapplicable in the immediate vicinity of the tran- 
sition and renders Eqs. ill-defined at some points. 
However, the overall conclusions about the switching 
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events will remain the same as long as the interval of 
inapplicability is small enough. 

We also find that the overdamped planar approxima- 
tion does not work well when a saddle point of magnetic 
energy is stabilized by spin-transfer torque, e.g. during 
the operation of a spin-flip transistor!^ Description of 
such cases in terms of effective planar equations requires 
additional investigations. 
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APPENDIX A: VECTOR DEFINITIONS 
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FIG. 3: Definitions of the tangent vectors and polar angles. 

We use the standard definitions of polar coordinates 
and tangent vectors (see Fig. [3]): 

n = (sin 6* cos 0, sin 6* sin cos 6*) 
ee = (cos0cos0, cos6'sin0, — sinf?) (Al) 
60 = (— sin 0, cos 0, 0) 

When 9 = TT /2 + 56 a linearization in 50 gives 

n « (cos sin 0, — (5(?) 
ee ~ -{S9cos(t),59sm<j),l) (A2) 
60 « (— sin (/), COS 0, 0) 

For two unit vectors n^'^\ i = 1,2 with polar angles 
{9i,(j)i) the scalar product expressions are 

(n'-'-' ■ Bg*^ ) = sin 9j cos 6i cos((/)j — (j^i) — cos 9j sin 9i 

(n(j") . ) = sin 9, sin(0, - 0,) (A3) 
Linearizing (jA3P with respect to small 59i for arbitrary 
values of 9j one gets: 

(n^^^ • Cg*') w — sm9j50i cos{(f>j — 0i) — cos^j 
(n(^) • ej*)) « sin^, sin((/), - (A4) 

Linearization of (IA3P with respect to both 59i and 69j 
gives 

(n(^') • 4'^) « -59, cos(0j - 00 + 59 j 
(n(^).e«) « sin(0,-0,) (A5) 
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